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Abstract 

An exact solution of the homogeneous wave equation, which was 
found previously, is treated from the point of view of continuous 
wavelet analysis (CWA). If time is a fixed parameter, the solution rep- 
resents a new multidimensional mother wavelet for the CWA. Both the 
wavelet and its Fourier transform are given by explicit formulas and 
are exponentially localized. The wavelet is directional. The widths of 
the wavelet and the uncertainty relation are investigated numerically. 
If a certain parameter is large, the wavelet behaves asymptotically as 
the Morlet wavelet. The solution is a new physical wavelet in the 
definition of Kaiser, it may be interpreted as a sum of two parts: an 
advanced and a retarded part, both being fields of a pulsed point 
source moving at a speed of wave propagation along a straight line in 
complex space-time. 



Keywords: physical wavelet, acoustic wavelet, localized wave, pulse, wave 
equation, mother wavelet, continuous wavelet analysis 
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1 Introduction 



The wavelet analysis has a great many different applications in signal and 
image processing (see [I], [2]), in physics and astronomy (see [I], [5], [6]). It is 
also used for developing efficient numerical algorithms for solving differential 
equations [7], [H]; however, mother wavelets are usually not associated with 
the solutions of differential equations under consideration. 

There is also an analytic approach to the problems of wave propagation 
proposed by Kaiser [9] , where the technique of wavelet analysis is developed 
for the decomposition of solutions of the wave equation in terms of localized 
solutions, which are called physical wavelets. They are constructed by means 
of a special technique of analytic continuation of fundamental solutions in 
complex space-time and can be split into two parts: an advanced fundamental 
solution and a retarded one. The physical wavelet as a localized solution 
of the wave equation has also been given in [10J. Applications of physical 
wavelets were discussed in [TTj-fTo]. 

In the present paper, we treat a new wavelet, which is at the same time 
a localized solution of the homogeneous wave equation in two or more di- 
mensions. This solution has been previously found and discussed in [IT], 
[T8] and generalized in [19]. It was named the Gaussian Wave Packet. We 
study its properties from two points of view. First, the solution can be taken 
as a mother wavelet for continuous wavelet analysis if time is a parameter 
and can be used in signal processing without being connected with any dif- 
ferential equation. Secondly, this solution should be regarded as a physical 
wavelet, i.e., it is an analytic continuation to the complex space-time of the 
sum of advanced and retarded parts of the field of a point source moving at 
a speed of wave propagation along a straight line and emitting a pulse that 
is localized in time. It is natural to decompose nonstationary wave fields in 
terms of these solutions, using the techniques of wavelet analysis. 

The aim of the paper is a detailed investigation of wavelet properties of 
the Gaussian Wave Packet for a fixed time and its properties as a solution 
of the wave equation in view of its further application to problems of wave 
propagation. For example, the decomposition of the solution of the initial 
value problem for the wave equation in terms of wavelets has been proposed 
by us in [20] . 

In Section [21 we give a brief review of the main facts of continuous wavelet 
analysis in one and two dimensions. 

In Section [31 we show that the Gaussian Wave Packet for a fixed time can 
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be regarded as a wavelet, give some estimates of it, and present its Fourier 
transform. We show that both the wavelet and its Fourier transform have 
an exponential decay at infinity. The wavelet has not only zero mean but all 
zero moments as well. 

In Section HI we discuss the asymptotic behavior of the Gaussian Wave 
Packet as some of the free parameters become large. We compare the packet 
with the nonstationary Gaussian Beam [21] , i.e., with the solution of the 
wave equation localized near the axis. We give the Gaussian asymptotic of 
it reducing it to the Morlet well-known wavelet [5j. 

In Section [5l we discuss the results of numerical calculations of the centers 
and widths of the packet in both the space and spatial frequency domains. We 
specify how fast these characteristics tend to asymptotic ones, with respect 
to the Morlet wavelet. We investigate the Heisenberg uncertainty relation 
for this wavelet, depending on the parameters and check how far from the 
saturation it is. We also obtain results for the nonasymptotic case where 
the wavelet corresponds to the solution that describes the propagation of the 
wave packet of one oscillation. This case may find applications in optics. We 
specify when the wavelet is directional [3J and calculate its scale and angular 
resolving powers. 

Section [6] gives a generalization of the above results to the case of an 
arbitrary number of spatial dimensions. 

In Section [71 we establish an analogy between the new wavelet and phys- 
ical wavelets of Kaiser [PJ. We show that it may split into incoming and 
outgoing parts, each solving a nonhomogeneous wave equation. As a source, 
we take new one-dimensional time-dependent wavelets, moving in the com- 
plex space-time at a speed of light. 

2 Main formulas of continuous wavelet 
analysis 

Wavelet analysis is a method for analyzing local spectral properties of func- 
tions (for example, see [I] - [1]). Wavelet analysis also allows one to represent 
any function of finite energy as the superposition of a family of functions 
called wavelets derived from one function called a mother wavelet by shifting 
and scaling its argument in the one-dimensional case and also by rotating it 
in the case of several spatial dimensions. By analogy, the Fourier transform 
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represents a signal as the superposition of oscillating exponents derived from 
one exponent by changing its frequency. 



2.1 One-dimensional wavelet analysis 

Let us give a brief review of some basic facts concerning the wavelet analysis 
of functions dependent on one variable x (for more detail, see [I] - jl], 0). 
Let a function (p(x) have a zero mean, and let it decrease as |x| tends to 
infinity so fast that tp(x) E Li (R) f] L2 (R) . It must oscillate to be nonzero 
and to have the zero mean. We call such a function a 'mother wavelet' 
because we derive a two-parametric family of functions from it, using two 
operations that shift the argument by b and scale it by a: 

^ a ' b (x) = j-^ip ^——j , b E (—00, 00), a E (—00, +00). (1) 

Thus any function (p a ' b (x) from this family again has the shape of (p, but 
shifted and dilated. By means of these operations, we can 'place' <p(x) at 
any point of the x axis and change its 'size' to any size by the parameter 
a. Then we define the wavelet transform W(a,b) of any signal f(x) by the 
formula 



W{a,b) = J dx' f(x')p a > b {x'), (2) 

R 

where the bar over (p denotes complex conjugation. 

One of the best-known mother wavelets is the Morlet wavelet, which is 

ip(x) = exp (^--^^j [exp(— z/cx) - exp(-K 2 a 2 /2)] . (3) 

It is the difference of a Gaussian function, filled with oscillations, and a term 
that provides the zero mean of (f(x) and that is negligible if kg ^> 1. It is clear 
from formulas (CQ) - ([3]) that |VT(a,6)| 2 , defined by (j2J), provides information 
about the frequency content of the signal / in the vicinity of the size aa 
of the point b, and k/cl plays the role of a spatial frequency. So we may 
regard the wavelet transform as a window transform, the size of a window 
changing for different frequencies. Changing the size of the window makes 
the wavelet transform more precise as compared to the Window Fourier (or 
Gabor) transform. 
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We can also reconstruct the signal f(x) from its wavelet transform W(a, b), 
or, in other words, represent the signal f(x) as a superposition of elementary- 
signals ip a ' b (x). Moreover, the mother wavelet used for the reconstruction 
of f(x) may differ from the one used for the analysis. The reconstruction 
formula looks like this 

oo oo 

m = 7^- \% [ dbW(a,b)P*{x), (4) 
J a z J 

-oo 

where ip(x) is another mother wavelet, and the constant reads 

oo — 

c ^ = j dk mm, ( 5) 

— oo 

where the symboPdenotes the Fourier transform. If we use the same wavelet 
for the transform and reconstruction, we should put ip(k) = <f>(k) and get 
|(J5(A;)| 2 in this formula to calculate the coefficient C w = C v . 



2.2 Wavelet analysis in two dimensions 

Wavelet analysis can also be defined for the case of more than one dimension 
(see [5], [H], [15]). A mother wavelet in the case of two dimensions r = (x, y) 
is a function ip(r) G Li(M 2 ) IL2 (M 2 ) that has zero mean. The Morlet wavelet 
in two dimensions reads 

~ 2^ J [exp(-iKx) - exp(-K 2 a 2 /2)] . (6) 

We define a family of wavelets from the mother wavelet, introducing rotations 
as well as dilations and the vector translations as follows: 



a,a, b{r) = _ M -l_ M -l = — — . (7) 

a \ a J \ sma cos a J 
The wavelet transform is defined as 

W(a,a,b) = J d 2 r' f{r') <p a ><*> b (r'), r' = (x',y'), d 2 r' = dx'dy'. (8) 

H 2 
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Then the reconstruction formula takes the form 

oo 2ir 



f{r) = (hJ S / ' d2b I daw ( a > a > b ^ a,a ' b ( r ^ ( 9 ) 



R2 o 

where 



C,. - / d 2 k ^ k ^ k \ k = (k x ,k y ), d 2 k = dk x dk 7 , (10) 



and the Fourier transform ip(k) is 

ip(k) — J d 2 r ip{r) exp {—ik ■ r), d 2 r = dxdy. (11) 



3 New two-dimensional wavelets 

We consider a family of functions ^(t*) of two spatial variables x, y containing 
arbitrary real parameters t, v and positive parameters p, e, 7 : 

,p(r) = J - ^d=p?=L, r = {x,y), (12) 



vr y^x + ct — is ' 
s = y/1 -i9/j, (13) 

y 2 

= x - ct + , (14) 

x + ct — IS 

where K v is the Bessel modified function (MacDonald's function) [23] . The 
branch of the square root in formula (fTBl) with positive real part is taken. 
The choice of the branch of the square root in the denominator of (fl2"|) is 
not important, for the sake of definiteness we assume that it has positive 
real part. We intend to show that each of the functions from the family (1120 
is suited for the role of a mother wavelet with good properties. The same 
is valid for their derivatives of any order with respect to spatial coordinates 
and time. 

Function f|T2l has appeared in [19] in connection with the linear wave 
equation 

iptt - c 2 (ip xx + ijjyy) = 0, c = const. (15) 
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If we regard the parameter t as time, formula (TT2"|) gives an exact solution 
of (TO)]) , which is well localized if p 1 (see below). If v = ~, formula (TT21 
yields 

exp(-ps) 

^( r ) = / ; i • ' ( 16 ) 
V ^ + ct — 2£ 



If in addition e = 7, it leads to the exact solution of (U5l) which was first 
reported in [T7] and discussed in detail in [IS] . 

In this section we view ^ as a two-dimensional mother wavelet with t 
being a parameter, ignoring that it is a solution of (fT5j) . According to [I], 
[5], this is possible, provided that the following conditions are satisfied: 

J d 2 r \ip{r)\ < 00, J d 2 r |?/>(r)| 2 < 00, d 2 r = dxdy, (17) 

R 2 R 2 

i.e., ^(r) G Li(M 2 ) f"|L 2 (]R 2 ), and it has zero mean, i.e., 

d 2 rij(r) = 0. (18) 



To prove (fTTI) we note that formulas (|T3|) . f|T4l) imply that Re(s 2 ) > 1 and thus 
we get Re(s) > 1, |arg(s)| < 7r/4. Therefore, (Tl2|) has neither singularities 
nor branch points for real x,y and t. It is a smooth function of x,y,t, and 
its derivatives of any order with respect to x,y,t are also smooth functions. 
We also obtain 

, 2|2 , (x-ct) 2 + y 2 y\( x -ct) 2 + y 2 + 2e 1 -e 2 -A{ct) 2 ) 

|S 1 7 2 7 2 [(x + ct) 2 + e 2 ] ' 1 1 

If x and y are large to an extent that the third term in (fl9l) is positive, then 

|* 2 | 2 > (X ~j 2 + ?/2 . (20) 

7 

Hence for large x and y the Bessel modified function can be replaced by its 
asymptotics, resulting in 



{vsY l l 2 exp (— ps) 

^ (r) = / t ; ■ — 

V + cr — 2£ 



1 + 



1 



(21) 
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Noting that Re(s) > \s\/y2 because | arg(s)| < 7r/4, we conclude that ip(r) 
has an exponential falloff and ip G L x f] L 2 . 

To check that the condition f|T8l) is satisfied, we calculate the Fourier 
transform of ip{r) 



d 2 r%jj(r) exp (— ifc • r), fc = (/c x , A; y 



(22) 



E 2 



The calculations yield (see Appendix 1) 



r> 2u 

il){k) = [k(k + K) v+l/2 ] 



-i 
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x exp 



,7 



-(k + k x )^--(k-k x ) 



— ikct 



(23) 



where k 



2 2 7 (/e + k K ) 
shows that ip(k)\k =0 = 0, owing the term 



|fe|. The formula 

in the exponent containing the denominator k + k x . Therefore, the function 
ip(r) has zero mean ffl8|) . Conditions (TTTT) . f|T8l) enable ip{r) to be a mother 
wavelet. 

Moreover, the following relation holds: 



Ql+r. 



(24) 



fc=0 



for any integer nonnegative j, /i, /, and m. This condition, the smoothness 
and the exponential falloff of ip{r) mean that any derivative of ip may be 
viewed as a mother wavelet and that all the moments of the wavelet ip and 
its derivatives vanish, i.e., 



0. 



03 x o^y 



0. 



(25) 



This property indicates that such wavelets could be useful in singular fields 

The wavelet ip has simple asymptotics for large values of p, which is dis- 
cussed in the next section. Below we present calculations of the wavelet ffl2l) 
and its Fourier transform (T23]) for moderate values of p when no asymptotics 
can be applied. It should be mentioned that the wavelet ffT2]) represents a 
wave of one oscillation when ^fp ps 1 (see Figures HJ [2]). This case is applicable 
in optics in the case of propagation of short pulses. 
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Re ip 



0.8^ 

A 




(a) Plot of the real part of tp(r). 
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0.4- 




(b) Plot of the Fourier transform V>(fc). 

Figure 1: An example of the new wavelet and its Fourier transform for p = 
0.5, e = 1, 7 = 0.25. 
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ReV> 




(a) Plot of the real part of ijj(r). 




(b) Plot of the Fourier transform V>(fc). 



Figure 2: An example of the new wavelet and its Fourier transform for p = 
1, e = 16, 7 = 0.5. 
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4 Simple asymptotics of the wavelet 



In this section, we study asymptotic properties of the Gaussian Wave Packet 
when the parameter p is large. First we note that we can replace the Mc- 
Donald function in (112j) by its exponential asymptotics (see [23]) and then 
we get (BID f° r an y x,y,t if p ^ 1. To prove this, we show that Re(s) > 1, 
then Re(ps) > p ^> 1, which provides that \ps\ ^> 1 and the exponential 
asymptotics is suitable. 

We intend to show first that the modulus of the exponent in fT2Tl) has 
a maximum at the point x = ct, y = 0. Let s = a + ib, then the relation 
2a 2 = (a 2 - b 2 ) + (a 2 + b 2 ) yields 

2[Re(s)] 2 = Re(s 2 ) + ^[Re(s 2 )} 2 + [Im(s 2 )] 2 . (26) 
Hence [Re(s)] 2 > Re(s 2 ), and taking into account that first 

2 

Re( s 2 ) = 1 + — — %- — (27) 

y ' 7 [ (x + ct) 2 + e 2 ] K ' 

and secondly 

\s 2 \ 2 = l+ (x ' 2 Ct)2 , (28) 

7 

which follows from (Tl9|) if y — 0, we get Re(s) > 1 outside the point a; = ct, 
y = 0. Formulas ffl3l . ffT4l) show that s — 1 if x = ct, y — 0, which is the 
point of the maximum of the modulus of exponent (1211) . Hence Re(s) > 1, 
next Re(ps) ^> 1, and the asymptotics fT2"Tl) can be used. 

Now we study the behavior of the Gaussian Wave Packet near its max- 
imum and show that it can be approximated by a nonstationary Gaussian 
beam. Consider a domain near the point x = ct, y = which is of order 

— 4 0(JLV !' of'V i< a <i. (29) 



7 VPV v^7 VP"/ 3 2 



We prove below that the wavelet ip(r) has a uniform asymptotics in this 
domain as p — > oo : 



ljj(r) = ^beam(K, r) exp 

47 



C [1 + O (p- 3a+1 )} , (30) 
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where k = p/(2j) and V'beam is a nonstationary Gaussian beam [21], [22] 



. . exp (% k o) . . 

^beam(^,r) = (31) 

Va; + cr — z£ 

C = p l/ - 1 / 2 exp(-j 5 ). (32) 

To prove (I3"0"j) . we decompose s from (I13p into powers of 9/j, which are 
of order 0{p~~ a ) by the formula (1251) and the inequality |y 2 /(x + ct — ie)| < 
y 2 /e. We obtain 

s = v^l - z0/ 7 = 1 - z^- + ^ + fi. (33) 

27 07^ 

We insert (133]) into (12T1) . substituting the expression for 9 ([Tlj) and the 
quadratic term # 2 /(87 2 ) = (x — ct) 2 /(87 2 ) + into it. It is easy to show 
that 

exp[-p(£i + £2)] = 1 + (p- 3a+1 ) . (34) 

We also note that the multiplier s u of the McDonald function in formula (fT21) 
has an estimate s v = 1 + 0(p _a ) and 0(p~ a ) = 0(p~ 3a+l ), because we have 
restricted a to the interval (|, |). The asymptotic formula (1301) is proved. 

Let us discuss formula (J3T?j) . The intervals where (x — ct)/j and y/^/ey 
vary are small compared to unit, according to ( 1291) . However they are large 
enough to ensure the exponential falloff of the function ip(r) on the edges of 
these intervals. To clarify this we note that the exponential terms in formulas 
( 130|) . (l3Tj) are of order 0(p~ 2a+1 ). They tend to infinity as p — > 00 by our 
choice of a < |. 

Asymptotics (1311 is expressed in terms of the field of the nonstationary 
Gaussian beam V'beam(«, v) [21]. If t is time, V'beam(^, v) is an exact solution 
of the wave equation (1T5|) with infinite energy. It is localized near the axis ox 
because Re(iK0) = —Key 2 /[{x + ct) 2 + e 2 ] and the cross section of the beam 
attains its minimum when x = —ct. Formula ( l30l) gives the Gaussian beam 
( 13T1) multiplied by the cutoff function exp[— k(x — ct) 2 /(47)]. The greater 7, 
the more elongated the essential support of ip, the closer it approximates the 
Gaussian beam (13T|) . 

Formula (130]) allows an additional simplification if 7 < e and 2ct/e = 
0(p~ a ). In view of estimates (1291) and the decomposition 

V ( V 2 \ V ( y 2 



7 \x + ct — ie J 7 \— ie [1 + i(x — ct)/e + 2ict/e] 
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iv t_ + o ( V 3Q+1 U o ( P - 3a+1 ) , 



57 Ve 

the asymptotic formula (!30|) takes the form 



^0) 



where 



C 



-is) 1 /' 



■ exp 



IK[X 



Ct) 



(x — ct) 2 y 2 



2ol 



2a 2 



[l + 0(p 



-3«+l 



(35) 



(36) 



4 7 2 /p, ^ = 1£ /p. (37) 

It is the Morlet wavelet (see pQ, [5]) with center at the point x = ct, y = 0. 
The numerical essential support of the wavelet, i.e., the points in the (x,y) 
plane where ip(r) is not negligible numerically, is an ellipse with semiaxes 
proportional to o~ x and cr y . The ratio <J x /a y 



2yl/ £ rules its asymmetry. 
We recall that, unlike the Morlet wavelet, the wavelet ip has zero mean and 
vanishing moments of any order. 

Another asymptotic formula can be obtained from the formula (!30l) if 
7 > e and ct/e > 1, but it is not necessarily small. This asymptotics was 
found and studied in [18] in detail with the restriction j — s. It is equal to 



ip(r, t) ~ exp 



IK\X 



Ct) 



[X 



ctf 



2a 2 



2aj 



+ i 



2ctK y 2 
AcH 2 + e< 



C 



(2ct 



where a 2 is the same as in (136 



ie)V2' 
(38) 

and is given in ( 1371) . The width of the packet 
along the axis oy is a function of time t and reads as a 2 = (T 2 {1 + Ac 2 t 2 /e 2 ). 

This asymptotics is obtained using fl33|) and decomposing l/(x + ct — ie) 
inside 9 in powers of x — ct. The quadratic terms in x — ct and y are taken 
into account. 



5 The uncertainty relation and directional 
properties 

In this section, we discuss numerical properties of the Gaussian Wave Packet, 
which are important for further applications of this new wavelet. We specify 
when the asymptotic case becomes valid. We also consider the wavelet in a 
nonasymptotic situation. 
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5.1 Widths of the wavelet ip and the uncertainty rela- 
tion 



Let us define the centers and widths of a function tp(r), r = (xi, ...,x n ). We 
denote by \\1pW2 the L 2 norm of ip: 

1/2 

d n r\ij(r)\ 2 \ , d n r = dx l dx 2 ...dx Q . (39) 



We define the centers and widths of a function as follows 



^=¥TiT2 I d n rx t \^r)\\ (10) 



1/2 

2 



Ax, = T - r / d n r (xt-xlY mr)\ 2 . (41) 



In the two-dimensional case we put n = 2, x\ = x and 22 = V- Using 
these formulas, we can also calculate the centers and widths of the Fourier 
transform of a function ip by replacing ip(r) by ip{k). We suppose that the 
essential numerical support of the wavelet, i.e., the points in the (x,y) plane 
where ip(r) is not negligible numerically, is the ellipse with semi- axes Ax and 
Ay. The essential numerical support of the Fourier transform ip(k) is the 
ellipse with semi-axes Ak x and Ak y and center at the point k = (k x , 0) (see 
Figure [3]) . 

The pairs of widths Ax, Afc x and Ay, Ak y satisfy the Heisenberg uncer- 
tainty relation, which holds for any function ip: 

1 A 7 A 1 

-, Ak y Ay>-. 

Equality holds only for the Morlet wavelet (see, for example [5]), which is 
the asymptotics of the Gaussian Wave Packet ffT21) as p — > 00. 

Our purpose here is to study the dependence of the widths and the means 
of the wavelet ip an d the uncertainty relation on the parameters. Taking e 
as the unit of distance, we rewrite the wavelet in terms of the dimensionless 
coordinates x' = x/e, y' = y/e, the time t' = t/e, and two parameters p and 
e/7 : 

#r0 = {ps) v K v {ps) 
\Jx' + ct' — I 



Ak x Ax>-, Ak y Ay > -. (42) 
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A k y 



2AA; X 



a 




2Ak x 



Figure 3: The essential numerical support of i/j. 
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x' - ct' + 



y 
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x' + Ct' 



nl/2 



(43) 



When p is large, the parameter p and the ratio e/ r y can be interpreted with the 
help of the widths of the wavelet Morlet flHBj) . which were denoted by <t x , cx y , 
and the mean of the longitudinal spatial frequency of its Fourier transform, 
which is k. The ratio e/7 = cr y /cr x characterizes the shape of the essential 
numerical support of the wavelet, the product k<j x = ^fp is the number of 
wavelengths on the width Ax, the product na y = a/ ke/2 is the number of 
wavelengths on the width Ay. 

We have calculated numerically the widths of the wavelet (Tl2|) , its Fourier 



transforms, and the left-hand side of the uncertainty relations (1421) as func- 
tions of p — 2/t7 in two cases. In the first case, the parameter k varies, 
the parameters e and 7 are kept constant. When p is large, it can be inter- 
preted as follows: the shape of the essential numerical support of the Morlet 
wavelet is kept constant and the number of wavelengths along the transverse 
and longitudinal widths increases with k. The results are plotted in Figures 
HI El In the second case, 7 varies, the parameters e and k are kept constant. 
For large p, this means that the number of wavelengths within the transverse 
width is kept fixed but the longitudinal width increases with increase of 7. 
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(a) Widths of the Gaussian Wave Packet tp(r) in relation to the 
widths of the Morlet wavelet, solid lines are for Ax/a x , dashed 
lines are for Ay/<r y : line (1) is for e/7 = 1/3, line (2) is for 
e/7 = 2/3, line (3) is for e/7 = 2. 




(b) Widths of the Fourier transform of the Gaussian Wave Packet 
ip(k) in relation to the widths of the Fourier transform of the 
Morlet wavelet, solid lines are for Afc x /(7k x , dashed lines are for 
Afc y /(7k y : line (1) is for e/-f = 1/3, line (2) is for e/7 = 2/3, line 
(3) is for e/7 = 2. 



Figure 4: The ratios of the widths of the wavelet ip to the width of the 
Morlet wavelet are plotted versus ^/p in spatial (a) and spatial frequency (b) 
domains for several values of e/7. 1R 
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(a) Widths of the Gaussian Wave Packet ip(r) in relation to the 
widths of the Morlet wavelet, solid lines for Ax/a x , dashed lines 
are for Ay/a y : line (1) is for 2ne = 4, line (2) is for 2ne = 8, line 
(3) is for 2ne = 64. 

1.8, 



1.4 



1 



0.6 



(b) Widths of the Fourier transform of the Gaussian Wave Packet 
ip(k) in relation to the widths of the Fourier transform of the 
Morlet wavelet, solid lines are for A/e x /(7k x , dashed lines are for 
Afc y /(7k y : line (1) is for 2ne = 4, line (2) is for 2ne = 8, line (3) 
is for 2ne = 64. 

Figure 5: The relative widths of the wavelet ip are plotted versus y/p in 
spatial (a) and spatial frequency (b) domains for several values of 2ks. 
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Figure 6: The left-hand side of the uncertainty relation for the Gaussian 
Wave Packet is plotted versus y/p, solid lines are for Ak x Ax, dashed lines 
are for Ak y Ay: line (1) is for e/'y = 1/3, line (2) is for 5/7 = 2/3, line (3) is 
for e/7 = 2. 
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Figure 7: The left-hand side of the uncertainty relation for the Gaussian 
Wave Packet is plotted versus solid lines are for Ak x Ax, dashed lines 
are for Ak y Ay: line (1) is for 2ks = 4, line (2) is for 2ks = 8, line (3) is for 
2k£ = 64. 
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The results are plotted in Figures El CD In all cases, the number of wave- 
lengths on the width Ax, i.e., y/p is laid off along the horizontal axis. It is 
the parameter y/p that controls the asymptotic behavior of the new wavelet, 
because the correcting term in fl36l) is of order 0(1/ y/p) when a = |. In 
view of comparing the widths with their asymptotics, the relatives widths, 
i.e., A x /cx x , and so on, are plotted as ordinates in Figures HI El 

Figure H] shows that all the widths tend to their asymptotic values as 
yjp — > oo, the parameter e/7 being fixed. The widths of the wavelet in the 
coordinate domain are larger than their asymptotics a x , a y (Figure Haj). The 
widths of the wavelet in the spatial frequency domain Ak x , Ak y may be both 
smaller or larger than the same widths of the Morlet wavelet o"k x = l/(2er x ), 
o"k y = l/(2er y ) (Figure l4bl . The larger the parameter e/7, the closer to 
the saturation the Heisenberg uncertainty relation (see Figure [6]) and the 
smaller Ay/cr y (Figure Haj) and Ak x /<Jk x ( Figure lib]) . However Ax/cr x and 
Ak y /<Tk y increase somewhat with increase in e/7, remaining smaller than 
their counterparts. 

To interpret Figures El [3, we note that the rate of convergence of the 
wavelet to the Morlet asymptotics is determined by terms of order p~ a, y /e = 
(2fi;) _a 7 1_Q /e, which must be small. However these terms increase with 7 if 
£ and k are fixed. 



5.2 Directional properties of the wavelet ip 

A wavelet is called directional (see [5] for greater detail) if the essential 
numerical support of its Fourier transform lies in a convex cone in k space 
with its vertex at the origin and the angle a at the vertex (See Figure [3]). 
Using the asymptotic formula fl5B|) . we can easily prove that the Gaussian 
Wave Packet is a directional wavelet when p is large enough. This follows 
from the fact that k x = k = pj (27), Ak x = yfpj (47) as p — > 00, 7 and e are 
fixed. Thus the inequality fc x > Ak x holds, which ensures that the origin lies 
outside the ellipse. 

We can calculate numerically the scale resolving power (SRP) and the 
angular resolving power (ARP) if the wavelet is directional. These quantities 
are especially important for numerical calculations, mainly in determining the 
minimal sampling grid for a lossless reconstruction of the image (for more 
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Figure 8: Scale and angle resolving powers for the Gaussian Wave Packet are 
plotted versus y/p, solid lines are for the SRP, dashed lines are for ARP: line 
(1) is for 5/7 = 1/3, line (2) is for 5/7 = 2. 











T) 
















\ 

























Figure 9: Scale and angle resolving powers for the Gaussian Wave Packet are 
plotted versus ^/p, solid lines are for the SRP, dashed lines are for ARP: line 
(1) is for 2ke = 4, line (2) is for Ike = 8, line (3) is for 2ks = 64. 
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detail, see [3J, The resolving powers are determined by the expressions 



SPRW = k i^k (44) 



'(k x ) - (Ak x ) 2 
APR{i[)) = 2 arccot v ^ = a. (45) 

As SRP tends to 1 and ARP tends to 0, the wavelet becomes more sensitive 
to small singularities and to angular details of the analyzed data [5]. The 
dependance of the angle and scale resolving powers was calculated by using 
MATLAB and is presented in Figure [8] and Figure [9] for the cases of varying 
k and 7, respectively 



6 Multidimensional case 

6.1 Definition and main properties 

Here we give a generalization of the wavelets constructed in previous sections 
to the case of many spatial dimensions. The counterpart of (fl~2l) in M. n is 



*(r) = M^M (46) 

^ y 1 {xi + ct-ie 2 y/ 2 (x l + ct-ie 3 y/ 2 - ...-(xi + ct-ien) 1 / 2 1 ; 

where r = (x\, x 2 , x 3 , . . . , x n ), p, 7, e 2 , . . . , e n are positive parameters, the 
function s is defined by (|T3|) . and it depends on r,t, and on the parameters 
via 9, which is as follows: 



e = Xl -ct+ ^ 2 — + ^— + ... + . (47) 

xi + ct — 1E2 x\ + ct — ie 3 x\ + ct — te n 



If t is time, the function ip(r) satisfies the wave equation 

ip tt - c 2 {ip XlXl +ip X2X2 + . . . + ip XnXn ) = 0, c = const. (48) 
If v — ~, formula (l4*6l) is transformed to an exact solution of ( HBl) . found in 



exp 

^ (r) = (27 + ct - ie 2 ) 1/2 (xi + ct - 2£ 3 ) 1/2 • . . . • (xi + c* - ^ n ) 1/2 ' (49) 
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The exact solution (1461) was mentioned in [19J. The Fourier transform of the 
wavelet (1461) is found in the Appendix. It reads 



2u 



x exp 



${k) = a n P — [k(k + hr+^ 2 ] 

k + h k\e 2 + kje 3 + . . . + kle n 
-7 



-i 



2k 2 7 



2 ' 2{k + k x ) k + h 

where k = (ki, k 2 , . . . , k n ) is a wave vector, k — \k\, 



ikct 



(27T 



,n/2_i7r(n-l)/4 



The axially symmetric wavelet follows from ()46l) . (1471) if e 2 = £3 
e n = e. It reads 

^ { P sfK v { P s) 



Its Fourier transform is 



7T (xi + ct- ie)^- 1 )! 2 ' 



(50) 



(51) 



(52) 



^k) = a n P — [kik + hY^-^ 2 ] L x 



2i> 



1 -1 



x exp 



(53) 



All the properties of (Tl2j) are extended to (l46|) with minor corrections, i.e., it 
is a smooth function of coordinates, its moments of any order are zero. The 
same is true for its derivatives of any order with respect to coordinates and 
t. The modulus of PST) has a maximum at the point X\ = ct, Xj — 0,j — 
2, . . . , n. It has asymptotics fl30l) as p — > 00, where 



■0beam( K ) r ) 



exp (? k 1 



(x x + ct- ie 2 ) 1 / 2 ...{xi + ct- ie n y/ 2 ' 
which is uniform in the domain 



xi — ct 

7 



O 



1W3 



o 



j — 2, 3, 



(54) 



(55) 



and I < a < |. If the parameters 2ct/ej, p~ ar y/ej, j = 2, .., n, are small, the 
asymptotics of the formula (I3"U|) with ^beam from (154")) represents the Morlet 
wavelet 



ip(r) ~ exp 



z/t(xi — ct) — 



(x\ — ct)' 
2a 2 



. 1.- ■-) 

2a 2 



2^ 



x 
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c 

X 7 „V \M2( „v M/2 ? „V M/2' ^ 56 ) 



(-i e2 )i/2(-ie 3 )i/2...(-i en )i/ 
where of = 47 2 /p, erf = 7£ 2 /p, ••• o"n = 7 £ n/P- 

6.2 Coefficient 

For the usage of the Gaussian Wave Packet as a mother wavelet, we need to 
calculate the coefficient (for example, see PQ, [5J) defined by the formula 

i" 

Although we cannot calculate this integral analytically in the general case, 
we can simplify this expression. It should be noticed that 

o 

(ic) n+1 f 1 
n — / dtt n exp(-ikct) = — -. 
n\ J y ' k n+1 

—ioo 

This equation makes it possible to treat the integrand in (1ST)) in the way 



$(k)\ 2 e -«(»-l)/4 7 M)/2( ic )n+l 



\k\ n (2ir) n / 2 p 11 - 1 n! 



y dtt n $ 2 (k,t). (58) 



Here we denote the Gaussian Wave Packet with parameters (7, e 2 , ...e n , 
by -0 and the packet with parameters (27, 2^2, 2e n , «, 2i/+(n — 1)/2) by r/> 2 - 
Then we can view the integral (1ST)) as a Fourier inverse transform calculated 
at the point r = 0. Substituting the expression on the right-hand side of 
(jSHJ) into (1ST)) and changing the order of integrals, we get 

_ i7r (n-l)/4 ^(n-l)/2 (j c \n+l °r 

C * = - {2K)n/i V 1 n! / ** B ^ '*)- (59) 

—ioo 

In the case where n = 3 and 6" = 7, we can explicitly calculate the 
coefficient from the formula (1ST)) if the coefficient 2v is a nonnegative 
integer. We calculate this integral in the spherical system of coordinates and 



23 



use the integral representation of McDonald's function K\{z) (see [23] )• We 
denote k = \k\. As a result, we obtain 



d 3 k 



ifci 3 



(27T) 



3^J1 / J3i. 

2 - / ifc| 5 (|fc| + k x y»+ 2 



r 



exp 



-27|fe| 



4k 2 7 



v y 7 4 ^ ^„ m! 



4i/+m— 1 



-4f+m— 5 



m=0 



7" 



— K m+3 (4/«7). 



(60) 



7 New wavelets and pulse complex sources 

In previous sections, we discussed properties of a Gaussian packet viewed at 
a fixed moment of time as a wavelet. Now we are going to review properties 
of a Gaussian packet as a physical wavelet. The term "physical wavelet" has 
been first suggested by G. Kaiser in [9]. These wavelets have been obtained 
as a sum of fields of two point sources with a special time dependence. 

First we consider two wave equations with moving 5-sources on the right- 
hand sides: 

Uu ~ c2 (*4 + O = ±0(9, t) 5 (x + ct)5(y), 
<p(q,t) = A^qexp(-2iqct), A = -e~ in/4 4:^c 2 . (61) 
The equation, say, for u + can be solved if we seek a solution in the form 

u + (x,y,t) = exp(iqa)g(j3,y), a = x — ct, (3 = x + ct. (62) 

This leads to the Schrodinger equation for the function g(P,y): 

&m + g yy = -^e- l ^ 4 ^6(f3)6(y)- (63) 

The solution of this equation is merely a well-known fundamental solution 
for the Schrodinger equation. So we get 



u^{r,t) = — exp 

y/X + Ct 



iq x — ct + 



ir 



x + ct 



(64) 
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where 9 is the Heaviside step function. We interpret it as a field of a source 
moving in the negative direction of the x - axis at a speed c, emitting an 
elementary pulse <fr(q,t). The solution does vanish in front of the moving 
source, i.e., for x < —ct, this is why it is called retarded. We denote the 
advanced solution, which does vanish behind the moving point, by u~ . It 
reads 



u 



(r,t) 



0(-x - ct) 



exp 



iq I x — ct + 



x + ct 



(65) 



\/x + ct 

It can be interpreted as a field absorbed by the source moving in the negative 
direction of the x - axis at a speed c and producing the current —(pit) in this 
source. 

We note that u = u + + u~ is a sourceless solution, i.e., it satisfies the 
homogeneous wave equation, because 9(x + ct) + 6(— x — ct) = 1. However, 
it has a singularity at the point x = —ct. In order to eliminate it, we shift 
x by —is/2, e > and t by —is /{2c) to the complex plane. This leads to a 
nonstationary Gaussian beam fl3Tl) : 



u(r,t) = ip heSkm (q,r) 



1 



\Jx + ct — is 



exp 



iq I x — ct + 



x + ct — ie 



. (66) 



We get a Gaussian beam, using the sum of fields of moving sources, one 
of which emits and the other absorbs a pulse <p(q,t) that is proportional to 
exp(— 2ictq). This result makes it possible to use an arbitrary function of 
time t instead of an exponent on the right-hand side of flBTl) . decomposing it 
into a Fourier integral. We apply 



$(t) = J dqF{q)4>{q,t) = B a v ~ 1/2 K u _ 1 / 2 (<j), 



(67) 



a 



2k 7 (1 + 2zct/ 7 ) 1/2 , B = -4cV i7r/4 — , 



J-(q) = aq v exp 



-7 | q + — 

Q 



(68) 



(69) 



instead of <J)(q, t) as input source functions in (IBT!) and obtain (fl2]) instead of 

V ; beam(g, »*). 
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8 Conclusions 

We suggest that we should consider an exponentially localized solution of the 
wave equation in several spatial dimensions, found earlier, from the point 
of view of continuous wavelet analysis. We investigate its properties as a 
mother wavelet and as a solution of the wave equation in view of its further 
application to the study of local properties and singularities of acoustic or 
optic fields. 

We show that, depending on the parameters, the solution represents a 
short pulse of one oscillation, or a wave packet with the Gaussian envelop 
filled with oscillations, or a nonstationary Gaussian beam multiplied by a 
cutoff function. 

The solution for a fixed time is a multidimensional wavelet with all zero 
moments. Its Fourier transform is calculated explicitly, and it is exponen- 
tially localized. The widths of the new wavelet in the position domain and 
in the spatial frequency domain, and the Heisenberg uncertainty relation are 
numerically investigated. 
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A Fourier transform of the Gaussian Wave 
Packet 

A.l Fourier transform of the two-dimensional packet 

The Fourier transform of the wavelet ffl2|) can be calculated using the Fourier 
decomposition of the Gaussian Wave Packet f|T2l in terms of the Gaussian 
beams ( |3TT) with the help of the formula 



oo 




(70) 
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To obtain this decomposition, we use the formula (see [2"3] ) 



- (ps) v K y (ps) 

7T 



dqJ-(q) exp( i q i 



(71) 



JF(g) = aq v exp 



-7 9 + 



(72) 



1 p 



2;/ 



""v^FC^)"' (73) 
where the left-hand side ( 17T1) depends on 6* via s from ( 1131) . We note that 
this relation is valid in the upper complex half-plane Im(#) > 0, where 9 
lies when x, y and time t are real. Dividing both sides of (ITTj) by the term 
\/x + ct — «£, we obtain ATP]) . So the Fourier transform of the Gaussian Wave 
Packet is determined by the formula 



ip(k) = / dgJF(g) ^ beam (fc, g). 



(74) 



The Fourier transform of the Gaussian beam (13T1) is found to be 

,2 



^beam(fc) 



, 2 exp (— ik ■ r) 
a r exp 



\/x + ct 



IE 



iq x — ct + 



y 



x + ct — is 



dx exp(—ik x x) exp[zg(a; — ct)} I y (q, x), 



(75) 



where 



I y (q,x) 



dy - : — — exp ( —ik y y + iq 



s/x + ct — is 



x + ct — is 



e iT/i exp 



K 

—i-r-(x + ct — is) 
Aq K 



(76) 
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Inserting (I75|) into (175|) . we obtain 

oo 

^beam(fc) = dx exp(-ik x x) exp 



zo(x — ct) — i~r-(x + ct — ie) 

4g V 



27r exp 



-M9 + 



4f/ 



ct - 



gf 
4g 



5 [ -kx + q--r-) e*^ 4 



k 2 
Aq 



(77) 



In view of further calculations of the integral with respect to q, we modify 
the 5 - function in (1771) . The roots of the argument of the 5 - function, i.e., 
the roots of the equation 



- K + q ~ / = 
Aq 



(78) 



are qi^ = (k x ± k)/2, k — \k\. The negative root does not contribute to the 
integral, because T{q, e) = if q < 0. Taking into account the relations 



5 ~K + q 



Aq] k 2 + Aq 



Aq 2 ( k x + k 

Hi o — 



k -\- k x 



5 , 

Aq q 2(k + k x ) ' 2 
we obtain 



Aq 2 



q k x + k 



i>be&m(k) = B(k)5 q 



k 2 + Aq 2 k 



k x -\- k 



2k 



(79) 
(80) 

(81) 

2{k + k x )J 2k " \jk + k x ^ 

After inserting (jS2J) into ([71]) . the integral disappears and the resulting ex- 
pression contains J-[(k + k x )/2], which is as follows: 



B(k) = 2tx exp I —ikct 



k 2 e \ k + k x i 



in/4 



2tt 



T 



k ~\~ k x 



w+l 



Taking into account 



(k + k x y+ i 



2(k + k x ) 



exp 



k + k x 2k 2 7 
2 /c H - /u x 



(83) 



(84) 
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we obtain 



r(k) 



,2v 



${k) = 2ixe mj ^ V — [k{k + k x y +1/2 ] 1 
x exp -ikct -(k- fc x )- - (fc + fc x ) - - - + fc 



(85) 



A. 2 Fourier transform of a multidimensional Gaussian 
Wave Packet 

In the multidimensional case, the Fourier transform of a Gaussian beam must 
be re-calculated. We put e 2 = e 3 = ... = e n = e from (jl6l) for simplicity. 
Instead of I y , the formula (175]) will contain the product J X2 ■ I X3 . . . J Xn , where 
n — 1 is the number of coordinates that are transverse to the direction of 
propagation. We obtain a formula for J Xk , k = 2, ...,n, by replacing y by Xk 
and a; by X\ in (1761) . The analog of (1771) will contain the sum fcf . = 
fc 2 — /c 2 instead of the term A; 2 = k 2 — k 2 , which does not yield any corrections 
and also the factor exp (in(n — l)/4) (n / q)^ 1 ^ 2 instead of such a factor for 
n — 2. Therefore, fl85|) is modified as follows: 

-i 



x exp 



—ikct — (k — ki)- — (k + fci)~ 



2 v ' 2 fc + fci 



(86) 



where 

A = (2vr) n/2 e i(n - lW4 - . (87) 
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